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Linear and nonlinear properties of convection in binary 
fluid layers heated from below are investigated, in particular 
for gas parameters. A Galerkin approximation for realistic 
boundary conditions that describes stationary and oscillatory 
convection in the form of straight parallel rolls is used to de- 
termine the influence of the Dufour effect on the bifurcation 
behaviour of convective flow intensity, vertical heat current, 
and concentration mixing. The Dufour-induced changes in 
the bifurcation topology and the existence regimes of station- 
ary and traveling wave convection are elucidated. To check 
the validity of the Galerkin results we compare with finite- 
difference numerical simulations of the full hydrodynamical 
field equations. Furthermore, we report on the scaling be- 
haviour of linear properties of the stationary instability. 
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I. INTRODUCTION 

Convection in binary fluid mixtures heated from be- 
low [1,2] is described by balance equations for mass, mo- 
mentum, heat, and concentration. The diffusive currents 
of heat and concentration that enter into the two latter 
balances are driven by generalized thermodynamic forces 
according to linear Onsager relations. They give rise to 
the Soret effect — temperature gradients change con- 
centration — and to the Dufour efffect — concentration 
gradients change temperature. In binary liquid mixtures 
such as alcohol-water [3-10] or 3 He- 4 He [11-13] the Du- 
four effect is negligible. Most of the research activity in 
the field of convection in binary fluid mixtures has been 
focussed on these binary liquid mixtures. 

However, in binary gas mixtures the Dufour effect is 
so large that it typically dominates the convective be- 
haviour whenever the magnitude of the Soret coupling 
strength, i.e., of the separation ratio ij) [2] is not negli- 
gible small. The importance of the Dufour effect in gas 
mixtures has two causes: (i) The Lewis number L = D/k, 
i.e., the ratio of concentration diffusion constant D and 
thermal diffusivity being of order 1 in gas mixtures is 
about 100 times larger than in liquid mixtures, (ii) The 
Dufour number Q measuring the contribution to the gen- 
eralized thermodynamic forces in the linear Onsager re- 
lations from gradients of the chemical potential that are 
caused by concentration gradients can be estimated [14] 
to be Q ~ 20-40 in gas mixtures. Now, the Dufour effect 
changes the (dimensionless) equation of motion 



{d t + u • V ) T = ( 1 + Q L ij) 2 ) V 2 T - Q L ij) V 2 C 

(1.1) 

of the temperature field T in two ways [14]. The "diag- 
onal" term QLip 2 V 2 T reflects an enhancement of tem- 
perature diffusion of relative size QLip 2 . The "offdiago- 
nal" contribution —QLip'V 2 C describes the direct effect 
of gradients in the concentration field C on the tempera- 
ture field. Both contributions to (1.1) are large when the 
size of the Soret coupling ij) is not too small. 

The influence of the Dufour effect on the onset be- 
haviour of convection in binary mixtures was determined 
within a linear analysis [14] of the convective perturba- 
tions of the quiescent conductive state. Here we first 
extend the exact analytical linear results of Lee, Lucas, 
and Tyler [15] for the stationary instability. Then, we 
mainly investigate various nonlinear convective proper- 
ties and how they are influenced by the Dufour effect — 
in particular for gas parameters. We mostly use an eight- 
mode Galerkin approximation to describe convection in 
the form of straight parallel rolls subject to realistic hor- 
izontal boundary conditions. 

Similar models for binary liquid mixtures have to cope 
with two difficulties: Boundary layer phenomena caused 
by the smallness of the Lewis number, L = O(10 -2 ), in 
liquids and the peculiar structure of the concentration 
field in traveling wave (TW) convection [16]. The re- 
stricted spatial resolution of a few-mode Galerkin trun- 
cation does not capture details of too fine a spatial field 
structure. Binary gas mixtures, on the other hand, are 
more favourable for such models: With L = 0(1) concen- 
tration boundary layer problems are less severe and the 
existence range of TW solutions is significantly reduced 
in parameter space since not only L but also the Prandtl 
number a is of order 1. In any case, we checked our 
analytical Galerkin results against finite-difference nu- 
merical test calculations of the full hydrodynamical field 
equations in order to assess the validity of the eight- 
mode Galerkin model. A positive feature of the latter 
is of course that it allows a convenient analysis of varia- 
tions with the control parameters Rayleigh number and 
separation ratio ij) and with the material parameters L, 
a,Q. 

The paper is organized as follows: In Sec. II the sys- 
tem, the equations, the boundary conditions, and the 
order parameters are described. Sec. Ill is preoccupied 
with analytical results of the stationary stability analysis 
of the quiescent heat conducting state. In Sec. IV we de- 
rive the Galerkin model and investigate the influence of 
the Dufour effect on linear properties of convective per- 
turbations and on the nonlinear solutions of stationary 
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and oscillatory convection. Sec.V contains comparisons 
of the model with linear and nonlinear results obtained 
from linear stability analyses and finite-difference numer- 
ical simulations of the full hydrodynamical field equa- 
tions. Sec. VI summarizes our results. In an appendix 
we present the corrected version of the linear stability 
analysis [14] for more idealized free-slip, impermeable 
boundary conditions. 



II. SYSTEM 



number L = — gives the ratio of time scales for concen- 
tration and heat diffusion. The Soret effect enters via the 
thermodiffusivity kx while 



1 dn(T,p,C) 

c p T 8C 



(2.5) 



quantifies the strength of the Dufour effect. In eq.(2.5) c p 
is the isobaric specific heat capacity and p, the chemical 
potential. 



We consider a fluid layer of height d between impervi- 
ous, perfectly heat conducting horizontal plates which is 
exposed to a homogeneous vertical gravitational acceler- 
ation g in z-direction. We impose a vertical temperature 
gradient so that the plates at z = are kept at tem- 
peratures To ± where To is the mean temperature of 
the fluid. The associated Rayleigh number 

agd 3 
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-AT 



(2.1) 



is given by the thermal diffusivity k, the kinematic vis- 
cosity v, and the thermal expansion coefficient 



1 8p(T,p, C) 

a = 

p dT 

The solutal expansion coefficient is given by 



/3: 



ldp(T,p, C) 

p 



dC 



(2.2) 



(2.3) 



where C denotes the concentration, p the pressure, and 
p the fluid's density. Throughout most of this paper we 
use dimensionless units which scale lengths by d, times 
by temperatures by and concentrations by . 



A. Equations 

The hydrodynamic field equations governing the sys- 
tem's dynamics are well known [14]. In Oberbeck- 
Boussinesq approximation they are 

V • u = (2.4a) 

(8 t + u • V ) u = - V (^- + g z) + v V 2 u 

Po 

+g [a (T - To) + (3 (C - C )] e z (2.4b) 

h 

To 

{d t + u • V ) T = k ( 1 + k 2 T L a ) V 2 T 
+k k T L a T V 2 C . 



(d t + u • V ) C = D V 2 C + D — V 2 T 



(2.4c) 
(2.4d) 



Here Co is the mean concentration of the mixture and 
D is the concentration diffusion coefficient. The Lewis 



B. Dimensionless deviations from the conductive 
state 

The stationary solution of the Oberbeck-Boussinesq 
equations (OBE) describing the state of pure heat con- 
duction without convection is 
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-*-cond — -^0 % 
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^cond — ^0 ~r ~ — Z 
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Pcond = p(z = 0) - pgz 



k T \ AT ' 
H T ) 2d 



(2.6a) 
(2.6b) 

(2.6c) 
(2.6d) 



Here Co is the mean concentration of the mixture. We 
pass over to reduced fields for the deviations from the 
conductive state 



(2.7a) 
(2.7b) 
(2.7c) 
(2.7d) 
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that obey the equations 

(d t -<tV 2 )v 2 w = <t (d 2 x + d 2 ) (e + c) + nlt 

(d t + u- V)0 = Rw+ (l + I/QV 2 ) V 2 

-LQtP V 2 c 
(<9 t + u-V)c = Ripw + IV 2 (c-ipe) 
V -u = 

NLT = e z • [V x V x ( (u • V ) u)] 



(2.8a) 

(2.8b) 
(2.8c) 
(2.8d) 
(2.8e) 



To derive (2.8a) we have applied twice the curl operator 
on eq.(2.4b). We have introduced the Prandtl number 
a = ^. The separation ratio 

V T a 

measures the Soret coupling. The Dufour number is 
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The Dufour effect enters into (2.8b) diagonally via the 
term LQip 2 V 2 9 that reflects an enhancement of heat dif- 
fusion — Q is positive — and offdiagonally via —LQipV 2 c 
which represents concentration-induced changes in the 
temperature field. Thus for small Soret coupling ip we 
can expect only small Dufour effects on linear and non- 
linear properties. 

The parameters L, a, ip, and Q depend on the mean 
temperature, concentration, and pressure of the fluid. 
Nevertheless, it is appropriate to characterize convective 
properties by L, a, ip, and Q instead of by the three ther- 
mal equilibrium quantities. Note, furthermore, that Q is 
known only poorly [17]. In order to select the range of 
parameters to be investigated here we used the follow- 
ing pieces of information: Hort et al. [14] have calculated 
Dufour numbers of order 10 using an ideal gas model, 
see also [17]. We limit ourselves mainly to the interval 
(0, 20). Since in gases concentration, heat, and momen- 
tum diffuse on the same time scales we mostly investigate 
mixtures with Lewis and Prandtl number 1. The sepa- 
ration ratio is considered in the interval (—1,0.25), that 
can be expected to contain the experimentally accessible 
range. 



C. Boundary conditions 

For a complete solution of the governing equations we 
need a set of boundary conditions for the three relevant 
fields w, 6, c. Except for the Appendix we impose re- 
alistic no slip, impermeable (NSI) boundary conditions. 
Since the concentration flux at the no slip boundaries is 
purely diffusive, 



-LV(C - ipT) for z = ±- 



we have to set 



d z (c - ip6) = forz = ± 7 



(2.9) 



(2.10) 



in order to avoid a vertical flux of solvent through the 
plates. In the conductive state the concentration flux 
vanishes identically. It is useful to introduce the field 

((x,y,z;t) = c(x,y,z;t) — ip 8(x, y, z; t) (2.11) 

instead of c(x, y, z; t) with the simpler boundary condi- 
tion 

d z ( = for z = ±^ . (2.12) 
The no slip boundary condition is described by 

w = = d z w for z = ±- . (2.13) 



Finally, since the temperature is fixed at the plates the 
deviation of the temperature from its conductive profile 
has to be zero for perfect conductors 



for z = ±- 



(2.14) 



More idealized, free slip, impermeable (FSI) boundary 
conditions are described in the Appendix. 

Using the (-field we get the following system of partial 
differential equations governing the convection in binary 
fluid mixtures 

{d t -aV 2 )V 2 w = a{8 2 x + 8 2 ) [(1 + V)0+C] +NET (2.15a) 
(d t +u- T7)6 = Rw + V 2 6-LQiPV 2 ( (2.15b) 
(d t + u- V)( = £V 2 C - tpV 2 8 (2.15c) 
V-u = (2.15d) 
NLT = e z • [V x V x ( (u • V ) u)] (2.15e) 

with the boundary conditions 

w = = d z w 



= d z ( for Z± = * 



(2.16) 



In eq.(2.15c) we have introduced an effective Lewis num- 
ber 



£ = L{1 + Qip 2 ) . 



(2.17) 



Therefore the Dufour effect is switched off by cancelling 
the term LQip'V 2 ( and replacing C by L. 



D. Order parameters 

To describe convection we shall use different order pa- 
rameters, (i) The maximal vertical flow velocity w max di- 
rectly measures the convective amplitude, (ii) The Nus- 
selt number 



TV" 



■^d z < e > x j 



= ±1/2 



(2.18) 



is the total vertical heat current through the layer re- 
duced by the conductive part R. Here the brackets imply 
a lateral average. To avoid the problem of determining 
the bulk heat current in the presence of a Dufour ef- 
fect we evaluate for convenience the vertical heat current 
through the fluid layer right at z = ±|: Not only heat 
advection but also any Dufour-induced contribution to 
the heat transport from vertical concentration currents 
vanishes at the NSI-plates. The reduced vertical heat 
current carried by convection alone, N — 1, measures the 
squared convective field amplitudes, (iii) Since w max and 
N — 1 do not characterize the concentration field we use 
the "mixing parameter" [16] 



M 



<{C- C ) 2 > 



x,y,z 



< (c c , 



Co) 2 > a 



(2.19) 
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M is the variance of the concentration field reduced by its 
value in the conductive state and thus characterizes the 
magnitude of concentration variations around the mean 
Co. In a perfectly mixed fluid, where all concentration 
deviations from Co vanish, M would be zero while in the 
conductive state with the Soret-induced concentration 
gradient M is defined to be 1. So 1 — M is an order 
parameter for convective states that is zero in the con- 
ductive state and approaches 1 for convection with per- 
fect mixing. Finally, (iv) propagating convection rolls are 
characterized by the oscillation frequency of the traveling 
convection wave. 



III. SCALING BEHAVIOUR OF STATIONARY 
STABILITY PROPERTIES 

As we will see later on the stationary stability thresh- 
old of the quiescent heat conducting state is the smallest 
one nearly all over the parameter space that is relevant 
for binary gas mixtures. Therefore, we compile, review, 
and extend in this section the exact analytical station- 
ary stability analysis of Lee, Lucas, and Tyler [15] that 
is based on the method of Chandrasekhar [18] (see also 
the paper of Gutkowicz-Krusin, Collins, and Ross [19], 
where, however, the diagonal contribution of the Dufour 
effect to the heat balance was ignored). Our analysis re- 
veals an interesting scaling behaviour of stationary stabil- 
ity properties that allows to scale away the Dufour effect. 
In addition to [15] we determine here the exact eigenfunc- 
tions. Furthermore, we present an analytical calculation 
of the zero wave number instability, an expansion of the 
critical Rayleigh number and the critical wave number 
around it, and further discussions, especially concerning 
the Dufour effect. 

Using 



(fe, 0). The linear equations of motion yield the relations 



R — RS ; S 
and 



£(1 + V)+V> 



/ ~ \ 1/3 

(i2/fe 4 ) (3.1) 



ik\fi 



(3.2a) 



1i 



k 

72 



V 1+T+T 2 +1+^ + i 



\/ l+T+T 2 -!- 



(3.2b) 



the stationary eigenfunctions on the marginal stability 
curve have the form 




cosh kz 
cos qoz 
cosh qiz 



cos kx ■ 



(3.3) 

Here the coordinate system is such that the horizon- 
tal wave number of the perturbation fields is (k x , k y ) = 



c 



R 



w 



-(i + V)0 
f c 

Lk 2 r y v> . 

R_ 1 -iy/3 f C 

Lk 2 r ' 2 V V> 



(3.4a) 
(3.4b) 

(3.4c) 



between amplitudes so that four unknowns 9, wo, wi, 
and u>l (6 and wq are real) and four boundary conditions 
remain. A solvability criterium yields the relation 

fetanh— |/ra |^\/3+ij gitanh^J + gotan^-| = 

p|gotan^- Im |^\/3 — gitanh^J — |gitanh^| 2 | (3.5) 

between r and k. It determines the marginal stability 
curve R s tab{k;p) depending on the parameter p given be- 
low. Then one can determine also the marginal ampli- 
tudes wo, wi, and w^ — with the normalization chosen 
such that 2(l+p)( = — T 2 k 2 p, they depend only on k and 
p. The nonlinear combination 



£(1 + V) £(1 + <?V> 2 )(1 + VO 



(3.6) 



of parameters L, Q, and ip that was called H by Lee, 
Lucas, and Tyler [15] is a scaling variable. Note that 
p vanishes in the pure-fluid-limit, ip = 0, of vanishing 
Soret effect. 

The marginal curve 



R s tab(k;p) = T 3 (k;p)k 4 



(3.7) 



obtained from solving (3.5) for r depends only via p on 
the parameters L, Q, and ip. Therefore, the critical wave 
number k c = k c (p) that solves dR sta b{k]p) / dk = is 
only a function of p. On the other hand, the critical 
Rayleigh number 



Rc = ^Rstab{k c {p);p) = ^R c {p) 
is a function of p and of 

£(i + y>) + y> _ v> i + p 

L ~ L p ■ 

The scaling relation (3.8) with £ = L(l + Qip 2 ) 
L ~ ( V 



(3.8) 



(3.9) 



(3.10) 



implies a significant simplification for practical calcula- 
tions of critical stationary properties: Only two func- 
tions, k c (p) and R c (p), have to be determined as func- 
tions of p to get k c and R c for all L, Q, and ip combina- 
tions. In Table I, we list, for several p, the scaled critical 
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Rayleigh numbers, R c (p), and the critical wave numbers 
k c (p). Using t c = (Rc/k 4 ) 1 / 3 (3.1), one then can deter- 
mine the critical vertical exponents go and gi (3.2). The 
vertical profiles of the critical eigenfunctions (3.3) can be 
obtained from wo and u>i in Table I. The amplitude £ 
appearing in (3.3, 3.4) is fixed by 2(l+p)£ = —T 2 k 2 p. 

Forp = 0, i.e., in the pure-fluid-limit ip = the scaling 
factor 5=1 and the scaling function has the critical 
value, _R°, of the one-component fluid: 

R c (p=0) = R° = 1707.762 and k c (p=0) = k° = 3.11633 

(3.11) 

The marginal stability curve shows the scaling behaviour 



R sta b(k; L, Q, ip) = —R sta b(k;p) 



(3.12) 



For values of p slightly below po 
we can calculate 



131/34 = 3.85294. 



3471468 p(^ -p) 



340023 p 2 + 2033552 p + 3779327 



(3.15) 



In Fig. 1 we show this k c (3.15) by a dashed line in com- 
parison to the exact result. Gutkowicz-Krusin et al. [19] 
and Knobloch and Moore [20] have given the expansion 
of Rgtabi'/L up to order k° thus obtaining the above 
stated value of 720. In (3.14) we present in addition 
the quadratic and quartic order and the p-dependence 
of k c whose determination uses the quartic order term in 
(3-14). 

In the literature dealing with binary liquid mixtures 
there are several expressions for the separation ratio ipo 
for which k c = 0. For zero Dufour effect, our calculation 
leads to 



where R s tab (3-7) is defined by the solution of (3.5). Thus, 
the Dufour effect can be scaled away in stationary stabil- 
ity proporties. Higher stability thresholds for odd eigen- 
functions can be obtained by replacing cos by sin and 
cosh by sinh in the terms of (3.3) containing the vertical 
spatial dependence. This results in an equation of the 
same form as (3.5), however, with tanh replaced by coth 
and tan by — cot. The case of heating from above with 
t < 1 and imaginary go can also be treated. Note that 
gotan(go/2) as well as gocot(go/2) are real for imaginary 

In Fig. 1 we show (a) the critical reduced wave number 
k c (p)/k® and (b) the reduced scaling function R c (p) / R° c , 
both versus p. For p — > —1, k c goes to about 7.48 and the 
scaling factor S in (3.8) goes to zero, which is the reason 
for the divergence of R c . At po = 131/34 ps 3.85 the 
critical wave number vanishes. We have evaluated po and 
the stability behaviour in the neighbourhood of po by a 
Taylor expansion of equation (3.5) up to order k 20 l 3 using 
the ansatz r = k~ 4 l 3 {a + bk 2 + ck 4 ) that follows from 
expanding the scaled marginal Rayleigh number R s tab as 

R stab (k) = T 3 k 4 = (a + bk 2 + ck 4 + 0(k 6 )) 3 

= a 3 + 3a 2 bk 2 + 3a(b 2 + ac)k 4 + 0(k 6 ) . (3.13) 

The value po for which the minimum of R s tab is located 
at k c = is determined by the requirement 3a 2 b = 0. In 
that case R s tab(k) increases proportional to k 4 . Perform- 
ing the expansion we get 



MQ = °) 



R ±-- 

^stab T — ^ , 

L 1 +p 



p n 

■H-stab 



2040 / 131 \ 

720 + 77p-( p -^r u 



1700115 
+ 1310309p 2 
+0(k 6 ) . 



, 2033552 3779327\ 4 



340023 



340023 J 

(3.14) 



f-L 



(3.16) 



with / = 1/po- Linz et al. [21] and Lhost et al. [22] have 
obtained from Galerkin approximations using free-slip, 
impermeable (FSI) and NSI boundary conditions values 
of / pa 1.62 and / ps 0.37, respectively. Knobloch and 
Moore [20] have extracted / ps 0.26 out of their numeri- 
cal stability analysis which has to be compared with our 
exact result of / = 34/131 = 0.25954.... 

Now we discuss the influence of the Dufour effect on 
the critical wave number's ^"dependence. Hort et al. 
[14] found, within their FSI Galerkin approximation, 
that with increasing Q the curve k c {ij)) formed a sad- 
dle at ij) = —1/3 when Q reached the value 27 for 
L = 1. This behaviour holds also for the exact re- 
sult with NSI boundary conditions. To show this, con- 
sider dk c /dij) = dp/dip ■ dk c /dp. Since fe c (p) (Fig. 1) is 
monotonous the extrema of k c {ij)) are given by the zeros 
oidp/dij), i.e., the roots of the third-order polynomial 



1 ,9 1 

-ib 2 = 

2 Y 2Q 



(3.17) 



One root is always greater than zero and not of interest. 
Two additional real roots first occur at Q = 27, ip = 
— 1/3 causing a saddle in k c {ij)) there. The appearence of 
this saddle is independent of L and a and holds for the 
exact NSI stability analysis as well as for the model [14]. 

We finally should like to mention that an application 
of the method presented here to the oscillatory stability 
analysis requires more numerical effort [20,23] because 
the vertical wave numbers q are not simple third-order 
roots but solutions of a fourth-order polynomial depend- 
ing also on frequency w. Therefore, exact analytical re- 
sults for the oscillatory threshold, the critical wave num- 
ber, and the Hopf frequency do not seem to be feasible. 
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IV. GALERKIN APPROXIMATION 



These modes drive nonlinear ones via the eqs (4.2): 



In this section we present our Galerkin model for re- 
alistic boundary conditions. Starting from this model 
we carry out a linear stability analysis of the conductive 
state before calculating nonlinear states, stationary as 
well as oscillatory ones. 



A. NSI mode truncation and model 

To describe convection in the form of straight rolls we 
truncate the spatial mode expansion appropriate to NSI 
conditions by 

w(x, z; t) = [wait) e~ ikx + c.c] Ci(z) (4.1a) 
6(x,z;t) = [9n{t) e~ ikx + c.c] V^costtz 

sin2-7rz (4. lb) 

((x,z;t)= [( 10 (t) e~ ikx + c.c.] + ( 01 (t)V2 cos wz (4.1c) 

where CAz) = co f^L - cos , Al , z 9 denotes the first 

^ / coshAi/2 eosAi/2 

even Chandrasekhar function with Ai = 4.73004. Since 
the modes #oi and (02 are linearly damped within an 
enlarged model's frame, as we have calculated before 
restricting the mode truncation to the above form, we 
do not display them explicitly. Furthermore, these two 
modes would violate a mirror-glide symmetry of the 6— 
and (-field that was found [24] for stationary and trav- 
eling roll patterns. The no-slip boundary condition is 
guaranteed by using the Chandrasekhar function. Im- 
permeability of the plates is ensured by confining the ( 
ansatz to modes with vanishing z-derivatives at z = ± | . 

Projecting the OBE (2.15) onto the eight modes con- 
tained in our truncation (4.1), we obtain the following 
generalized Lorenz model 



tX 



-5/X + 5ff[(l + V)Y+ ai U] 



(4.2a) 



tY = —q 2 Y + (r-Z)X + (3^k 2 V (4.2b) 
tU = -/3i£fe 2 U + g 2 V>Y + VX (4.2c) 



tZ 



X Y 



4^1 



tV 



with 



[auX-u + m^z + cv] 



;LQip , £ = L(1 + QV 2 ) 



(4.2d) 
(4.2e) 

(4.3) 



We used for the critical modes the following vector nota- 
tion 

X = {X 1 ,X 2 ) = ^^(Rew^Imw^) 

Y = (n,y 2 ) = 2a 3 &{Re6 11 ,Im6 11 ) (4.4) 

U = {U U U 2 ) = a 3 ^p(Re( 10 ,Im( 10 ) . 



2V2a 3 ^9 02 



and V 



C01 



(4.5) 



The constants are 

ai = 2ttA 2 {j^i ~ AJ-8I71- 4 ) 



a-2 
a 3 

ai 

Pi 
Vi 
Mi 



2Ai tanh 



4ttA 2 



^-A 2 tanh 2 ^ 



—V tanh ^ 

fc 02 +7T 2 
4 

3a s 
~2al 



0.4058 

0.6974 
0.5818 
-12.3026 
0.7585 

0.4930 
2.2916 
1.4770 



(4.6) 



The quantities k° = 3.098, q° = k° 
T = 4i = 0.05138 and b = ■. 



-7T 2 , R° c = 1728.38, 
2.0282 are critical 

properties of the model for ij) = [25]. a = 1.943c 
denotes a rescaled Prandtl number caused by the no-slip 
mode truncation [25]. In addition we use the reduced 
Rayleigh number r and the wave numbers k and q defined 

by 



_ _R i_k_ .2_r_ 

" ~ R° c ' *° ' q ~ qf 

as well as the quantities 

Af + k 4 -2k 2 a 4 



kf 



U02 



(4.7a) 



/ 



04 



a 4 Af + fc° 4 - 2k° c 2 a 4 



U02 



and 



04 



04 



(4.7b) 



so that k = q = f = g= l for k = This nota- 
tion is also used by Lhost et al. [22] who have examined 
the onset of convection in binary mixtures neglecting the 
Dufour effect and by Niederlander et al. [25] who have de- 
rived the analagous NSI model for one component fluids. 
Note, however, that in [25] the modes are reduced by k- 
dependent quantities — see, e.g., eqs. (2.8-2.10) in Ref. 
[25]. On the other hand, here we reduce the modes by 
critical properties of the ij) = reference system. Thus, 
here all fe-dependence of the system is displayed explicitly 
in the model equations (4.2). 

The Dufour effect influences the mode balances in two 
distinct ways: It causes a driving of the temperature 
modes Y and Z by the (-field modes U and V, respec- 
tively. The associated coupling strength ip = -\LQip 
vanishes when Q = 0. Furthermore, in the (-field equa- 
tions for U and V, the Lewis number L is replaced by an 
effective one, C = L(l + Qip 2 ). 

To facilitate the quantitative comparison with exper- 
iments or numerical solutions of the full field equations 
we evaluate among others the order parameters defined 
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in Sec. II D. Thus the maximal vertical flow velocity is 
given within the model by 



Ci(0) I X I = 12.20 | X 



(4.8) 



in terms of the amplitude of mode X. The reduced ver- 
tical convective heat current, evaluated at the plates, is 



N 



27iV2, 
R 



702 



z 

a 3 r 



(4.9) 



Thus the Nusselt number is related to the mode Z. As an 
aside we mention here a deficiency of the no-slip Galerkin 
approximation that was discussed in more detail in [25]: 
Since the velocity field is expanded in Chandrasekhar 
functions, i.e., a nontrigonometric basis the stationary 
vertical heat current is not z -independent. 

Into the mixing parameter M (2.19) the temperature 
modes as well as the (-field modes enter 



M 



1+ 



24 



r 2 V al^qf 
24 4 



r 2 tp 2 7r 2 
6 



U 2 + ■ 
64a 2 . ' 



-V> 2 Y 2 + 2V>Y • U 
, 2V> 



3a 3 



-ZV 



(4.10) 



B. Linear stability analysis 

We start the discussion of our model with the investi- 
gation of the linear stability of the conductive fixed point 
where all mode amplitudes vanish. The nonlinear modes 
Z and V do not couple linearly into the equations for X, 
Y, U and are damped away. Therefore we have to seek 
the stability thresholds of the matrix system: 



-5/ ag{\ + V>) oqax 
r -q 2 /3iV^ 2 

q 2 ip -PiCk 2 




(4.11) 

We calculate the stationary stability curve r stat (k) to be 

8QV> 2 

(4.12) 



r sta t(k) 



fg 2 tt 2 (1 + Qj, 2 ) 



1 + V 1 + 



«l 



The critical stationary wave number follows from 

= Ktat 6 + h Ktat 4 + h Ktat + fo (4-13) 

with the coefficients 



71-73 , 3V> 71 a i 



2C 



i> 71(71-73)"! 



fo 



7i72 , ip 7i "1(72 - 737i) 



2 2£ 



(4.14a) 
(4.14b) 
(4.14c) 



and the numbers 

71 = ^2 =1-029 ; 72 = -^ = 5.437 ; 73 = 2-^ = -2.564 . 

K c K c K c 

(4.14d) 

The critical stationary wave number as a root of (4.13) 
is a function of the scaling parameter p = only — 

different ip, L, Q combinations for which p is the same 
yield the same k c stat (p). Such a scaling behaviour was 
found analytically in the exact stationary stability anal- 
ysis of the full field equations by [15] and in Sec. III. How- 
ever, the marginal curve r stat (k) (4.12) of the model does 
not show the full scaling behaviour seen in Sec. III. 
The oscillatory stability curve is given by 



M[M) H) 



tptp 



{\ + i]>){\ + a) - a-Lip 



(4.15) 



and the Hopf frequency at r osc (k) 



4 



aiV> 



(1 + £)(? + £) 
(1 + + a) - 
(1 + V>)(£ -a)+ ai ij) 



(4.16) 



Here we have introduced wave number dependent fluid 
parameters 

a=l-a , £ = (3 1 ^-£ , and J = /3 * ^ . (4.17) 

Since the square of the critical oscillatory wave number 
is given by a root of a polynomial of at least degree ten 
we evaluated k c osc numerically by minimizing r osc (k). 

Fig. 2 shows the main results of our stability analy- 
sis: The Dufour effect destabilizes (stabilizes) the con- 
ductive state against the growth of stationary (oscilla- 
tory) convection and it shifts the critical curves r c stat {ij)) 
and Tg SC (ip) towards smaller ip. The critical stationary 
wave number k c stat {ij)) forms a saddle for Q = 27 and 
tp = — |. The wave number k c osc (ip) of critical oscilla- 
tory patterns decreases and the difference, k c stat — k c osc , 
increases with increasing Dufour effect. The Hopf fre- 
quency decreases with growing Q. An important fact is 
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that the V'-range of oscillatory instability sharply shrinks 
with increasing Dufour number. It remains an experi- 
mental challenge to prepare mixtures that have the right 
Q—ip parameter combinations to see this behaviour. All 
the above described properties are the same as those ob- 
tained from the full field equations with a numerically 
performed shooting analysis (compare, e.g., Fig. 7 of Ref. 
[14] with our Fig. 2). We refer to Ref. [14] for a more de- 
tailed discussion of linear properties which are not the 
main topic of this paper. 

Finally, we mention that the analytic oscillatory sta- 
bility analysis presented in [14] for idealized FSI bound- 
ary conditions contains a mistake. In the Appendix we 
present the correct formulae and a figure showing these 
results. We find that our corrected FSI results are closer 
to the exact NSI curves than the FSI results of Ref. [14]. 



C. Nonlinear convective states 

Here we elucidate the influence of the Dufour effect on 
porperties such as strength of convection, bifurcation be- 
haviour, heat flux, and concentration mixing in nonlinear 
states of stationary and oscillatory convection. 



1. Stationary convection 

The stationary solutions of our model representing 
steady overturning convection (SOC) in the form of 
straight rolls are given by 



X 2 
Y 



2 

-F 



mi V 4 mi 

r2i2 i Vl ~2 / 7 i -v-2 

— £ k + - — q ipip + X 
Mi 4/ii 



Z = —F 

U = ^*> 
Mi 



m «fv>V> + x 2 



4M 



V = Fiji 



Mi 



(4.18a) 
(4.18b) 
(4.18c) 
(4.18d) 
(4.18e) 



where we have introduced the following abbreviations: 



F ~ 1 = l ( " - ) M 2 ^ 1 + W + 

f \ 4 Mi Mi/ v ' 



4x 2 



i + V> 1 



Mi 



(4.18f) 



Mi 



1 + V> 1 



Mi 
_ aim 
Mi 



a = q H £ k -\ ipip(k + ——) 

4/3i 



(4.18g) 



(3 




T stat J 

(4.18h) 



The formulae (4.18) for the SOC solution are structurally 
similar to those of the analogous models [26,21] without 
Dufour effect and derived for idealized FSP [26] or FSI 
[21] boundary conditions. Also here, like in [26,21], the 
SOC solution does not depend on the Prandtl number. 

A positive X 2 bifurcates according to (4.18a) out of 
the conductive solution, X = 0, at the stationary thresh- 
old r st at where (3 (4.18h) goes to zero. The influence of 
the Dufour effect on the convective intensity can be seen 
in Fig. 3. There we display X 2 , the square of the verti- 
cal velocity mode, versus r for ip = —0.25, L = 1, and 
k = for several values of the Dufour number Q. For 
a fixed Rayleigh number r > 2 the strength of convec- 
tion is reduced with increasing Q. However, near onset 
it is enhanced: The Dufour-induced destabilization of 
the conductive state shifts the convective onset to lower 
values of r. The reduction of X 2 at large r is nearly pro- 
portional to Q. This can be checked by an expansion for 
r — > oo, which already holds for r > 3, where all curves 
X 2 (r) tend asymptotically to straight lines. 

The most conspicuous change in the SOC bifurcation 
behaviour with increasing Dufour effect is the gradual 
change from a strongly backwards bifurcation (Q = 
in Fig. 3) via a tricritical one to a forwards bifurcation. 
This behaviour is documented in a more global manner 
in Fig. 4: The bifurcation of SOCs with k = 1 is for- 
wards (backwards) in the shaded (white) L-i/'-region to 
the right (left) of the thick full curve of tricritical bifurca- 
tions. The shaded region strongly grows with increasing 
Q on cost of the white one. At the thin solid line of Fig. 4, 
the bifurcation threshold has moved to r stat = oo. So, to 
the left of it, the lower branch of X 2 (r) is disconnected 
from the X = solution. In this regime, convection 
branches out of the conductive state at a finite negative 
fstat, i- e., for heating from above. 

Note that the tricritical bifurcation line at, e.g., Q = 15 
in Fig. 4, is non monotonous — the Dufour-induced ap- 
pearence of strongly nonlinear L-ij) variations in convec- 
tive and in stability [14] properties is not surprising in 
view of the fact that the Dufour effect enters via LQij) and 
LQip 2 into the field equations. Thus, when decreasing 
here ip from zero, e.g., along the dotted line in Fig. 4, one 
can observe the succesion f—ft—fb—ft—ff—ft—fb 
of forwards (/), tricritical (t), and backwards (6) bifur- 
cations. 

In Fig. 3 we have seen that the flow intensity X 2 de- 
creases with increasing Dufour effect Q or with decreas- 
ing heating rate, r. In Fig. 5 we show that the structural 
changes resulting from either increasing Q or decreasing 
r are in fact the same: The simultaneous agreement in 
the structure of all fields — temperature, velocity, and 
concentration — is almost quantitative for the two com- 
binations Q = 10, r = 2.5 and Q = 0, r = 2 shown in 
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Fig. 5 and could be made perfect by judiciously choosing 
r-Q-combinations. Compare also the lateral profiles of 
these fields at midheight of the fluid layer as shown in 
the second row of Fig. 5. Thus, increasing both, Q and 
r, appropriately does not change the SOC state. 

In the second row of Fig. 6 we show, for various Soret 
coupling strengths ip, the influence of the Dufour effect on 
the r-variation of the model's Nusselt number (2.18,4.9) 
in SOC states. The influence of the Dufour effect on the 
Nusselt number is similar to that one on the flow inten- 
sity — cf. Fig. 3. First of all, the onset of convection is 
shifted with increasing Q to a smaller Rayleigh number. 
Simultaneously, at larger r, say above r = 2, the verti- 
cal heat current decrases with increasing Q. Both effects 
combined flatten the bifurcation curve of N — 1 versus r 
and/or change the backwards bifurcation topology at suf- 
ficiently negative ip into a forwards one. For example, in 
the extreme case of ip = —0.5, where for Q = the lower 
bifurcation branch (cf. dots in Fig. 6) is disconnected 
from the conductive state, already a Dufour effect of size 
Q = 5 brings down the onset r stat to about 3.47. Increas- 
ing Q further the bifurcation becomes forwards even for 
this strong Soret coupling ip = —0.5. 

In the third row of Fig. 6, the graphs of M (2.19,4.10) 
versus r show how the Dufour effect influences the mean 
square variation of the concentration field in the SOC 
states of the model. Remember that M is defined to be 
1 in the conductive state. Furthermore, the better the 
convective mixing of the fluid the smaller are concentra- 
tion variations and with it M . So we see in Fig. 6 that the 
Dufour-induced reduction of the convective flow intensity 
and of the Nusselt number at larger r is accompanied by 
a reduction of the convective mixing: For larger r, the 
parameter M increases with increasing Q. Roughly and 
qualitatively speaking the bifurcation behaviour of the 
flow intensity, w^ ax , an< ^ of the convective heat current, 
TV" — 1, are similar to 1 — M which measures the degree 
of convective mixing. 

2. Traveling wave convection 

Our Galerkin model has nonlinear convective solutions 
in the form of harmonic waves of constant amplitude 
traveling with constant phase speed v p = u>/k either to 
the left or to the right. The complex Galerkin modes 
(4.4) for this TW solution have the form 

X(<) = | X | e iut (4.19a) 
Y{t) = | Y | e i{ut+a ^ (4.19b) 
U(<) = | U | e^'+W (4.19c) 

with constant moduli | X |, | Y |, and | U | and 
phase differences a, (3. The real modes Z and V (4.5) 
are time independent as well. The hydrodynamic fields 
follow from (4.1). 



From the mode equations (4.2) one finds that both 
| X | 2 and the squared frequency w 2 (k) of the TW solu- 
tion vary linearly with the distance from the oscillatory 
threshold at r osc (k): 

|X| 2 =s TW {r-r osc ) (4.20a) 
w 2 = e>% + hw {r - r osc ) . (4.20b) 

The threshold values r osc (k) and Wjj(fc) as well as the 
slopes stw and fxw of the bifurcating TW solution 
branch depend on L, Q, ip, and a but the r dependence 
is always linear, i.e., of the same form found for Q = 
and liquid mixture parameters with FSI boundary condi- 
tions [21,27]. The subcritically (supercritically) bifurcat- 
ing TWs with stw < (stw > 0) are unstable (stable, 
at least close to onset). At 

r* = r osc - ^- (4.21) 

STW 

the TW solution ends by merging with zero frequency 
into the SOC solution (4.18). 

Fig. 7 shows existence boundaries of the TW solutions 
(stable or unstable). For a = 1 the model has TW solu- 
tions with k = in the Q-i/'-region below the thick lines 
for various Lewis numbers as indicated by the line type. 
These existence boundaries are determined by the merg- 
ing of the oscillatory and stationary bifurcation thresh- 
old, r osc (fe°) = r sta t(kc), for fixed wave number i.e., 
by the codimension-two (CT) condition with vanishing 
Hopf frequency wjj(fe°) = 0. Note that the existence 
range of TWs with the critical oscillatory wave number 
k c osc is somewhat wider than the one for TWs with as 
can be inferred also from Fig. 2. The boundary curves of 
Fig. 7 show that the ^-range in which TWs can be found 
is shifted to more negative values when (i) the Lewis num- 
ber increases, i.e., when concentration diffusion becomes 
more efficient or (ii) when the Dufour coupling increases. 
For L < 0.6 the grange of TW existence range splits 
into two pieces for sufficiently large Q. A similar feature 
was found in the linear analysis of the regions with non 
vanishing Hopf frequency (Fig. 6 of [14]). 

The TWs bifurcate supercritically (subcritically) out 
of the conductive state for the parameters in the shaded 
(white) regions of Fig. 7 below the respective thick 
curves. For parameters on the thin, vertically oriented 
lines, the TW bifurcation is tricritical. Obviously, for- 
wards bifurcating TWs disappear pretty soon with grow- 
ing Dufour number while backwards bifurcating ones can 
be seen for larger Q, provided the Soret coupling is suf- 
ficiently negative. The width of the shaded existence 
range in Fig. 7 of supercritical TWs, i.e., the distance 
I V'ct — ipt I between CT and tricritical separation ra- 
tio measured, e.g., for Q = and a = 1 is largest at 
L = 0.736. This distance decreases to zero for L — > 0. 

Our model also shows for positive ip TW solutions like 
the FSI model for Q = [21,27]. These TWs at ip > 
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branch in a secondary bifurcation at r* out of the SOC 
state whereas for ip < TWs bifurcate at the oscilla- 
tory threshold r osc out of the conductive state and end 
at r* in the SOC state. We should like to mention that 
numerical simulations of liquid mixtures gave no indica- 
tion for the existence of TWs at positive ip (cf. [16] for 
a short discussion) which suggests that they result from 
the mode truncation of the model. For gas mixtures nu- 
merical results are not available. Furthermore, the model 
TW states at positive ip appear in gases at large r (e.g., 
r > 10 for Q = 10, ip = L = a= k = \) where with 
the mode amplitudes being large the model's truncation 
approximation is presumably not justified. And, finally, 
these states occur for relatively large positive separation 
ratios that might be inaccessible experimentally. 

V. COMPARISON WITH RESULTS FROM THE 
FULL FIELD EQUATIONS 

In this section we compare linear as well as nonlin- 
ear properties of our eight-mode Galerkin approxima- 
tion with results obtained from the full hydrodynamic 
field equations. In Subsection V A critical linear model 
properties are compared with an analytical exact analy- 
sis of the stationary instability of Sec. Ill and with results 
obtained numerically with a shooting method for the os- 
cillatory instability [14]. In Subsection V B we compare 
nonlinear convective properties of the model with results 
from solving the field equations numerically with a finite- 
differences method. 



A. Linear properties 

First we compare the numerical exact critical wave 
numbers, stability thresholds, and Hopf frequencies of 
Fig. 7 of [14] with Fig. 2 showing in an analogous way 
the results of our model. Our approximation of the criti- 
cal wave number of the stationary patterns is nearly per- 
fect with errors less than 1 % in the whole Q— and ip- 
range investigated here. The stationary stability thresh- 
old r c stat itself is approximated with the same quality over 
the whole ^"interval for Q < 10 and for all investigated 
values of Q for ip > —0.3. The build up of a local min- 
imum and maximum in r c stat (ip) occuring for Q > 10 in 
our model starts for smaller values of Q in the numerical 
threshold [14]. This reflects small errors for moderate Q 
and strong ip < —0.6. Our model reproduces the V'-value 
where r c stat (ip) diverges with an accuracy of about 5 % 
since the model's stability curve does not show the the 
proper scaling behaviour (see Sec. Ill) resulting from the 
full field equations. 

The errors of our model for the oscillatory stability 
analysis are largely caused by a shifted ^"dependence. 
This shift from the model values to those of the shoot- 
ing analysis increases with stronger Dufour effect from 



for Q = to -0.3 for Q = 20 in the direction of the 
negative V^axis- After this transformation the model 
describes r c osc as well as the Hopf frequency with small 
errors. Even the slight kink in the curve Wjj(V') for Q = 5 
in the neighbourhood of ip ~ —0.5 is reproduced. In ad- 
dition, our model provides good approximations for k c osc 
if we restrict ourselves to values of Q less than 5. For 
stronger Dufour effect it cannot produce the strong ip- 
dependence seen in the shooting analysis, but it shows 
the main two qualitative effects of increasing Q on the 
critical oscillatory wave numbers: a decreasing value of 
k c osc and a building up of a minimum in the curve k c osc (ip). 

To conclude, our NSI-Galerkin-approximation is very 
good as far as the stationary stability analysis is con- 
cerned. The errors in its oscillatory part are small and 
increase with Q. But they are less relevant, say, from an 
experimental point of view because they occur mostly in 
a parameter range where the first instability is stationary. 

To get further insight into the reasons of the errors of 
our model we compare in Fig. 8 the vertical profiles of 
the critical modes for the stationary and the oscillatory 
instability, respectively, for Q = 5 as a representative 
example. Exact results are shown by full lines and the 
model's results by dashed lines. Because of the mirror 
symmetry of the critical modes at the mid plane z = 0, 
the modulus and phase profiles are displayed only in one 
half of the layer. 

Consider first the profiles of the stationary critical 
modes (left row in Fig. 8). Those of the #-field contain 
with increasing Q and \ip\ a small admixture of higher 
vertical modes beyond the cos -kz of the model. The crit- 
ical velocity field is fitted very well by the first Chan- 
drasekhar function excepted for a small error in its cen- 
tral maximum which increases with Q. The most impor- 
tant deviation occurs for the critical (-field: The approx- 
imation £ = const ignores the z-variation of the exact 
profile. But the constant lies for relevant fluid parame- 
ters within the range of variation of the exact result. 

In the right row of Fig. 8 we compare moduli and 
phases of the critical oscillatory modes of the model with 
the corresponding exact results. Also here, like for the 
stationary instability, the vertical profiles of the critical 
eigenfunctions of velocity and temperature resulting from 
the model agree quite well with the exact ones in modu- 
lus and phase. However, the model phases are constant 
and typically larger than the z-dependent exact phases. 
Again the largest deviations, in modulus and phase, oc- 
cur for the critical (-field. 

B. Nonlinear Properties 

To test the quality of the model's predictions for 
the nonlinear convective states we have performed some 
selected numerical simulations of the full hydrody- 
namic field equations (2.4) with a MAC/SOLA finite- 
differences method. This code has been employed suc- 
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cessfully for binary liquid mixtures [16]. 

1. Traveling wave convection 

For gas mixtures L = 1 = a and < Q < 20 we found 
in the full-field simulations no stable, large-amplitude 
TW solutions that have bifurcated subcritically in agree- 
ment with the model. Note that the existence range of 
nonlinear stable TWs on the upper TW solution branch 
in binary liquid mixtures at Q = rapidly shrinks to 
zero as L approaches 1 from below (Fig. 14 in [16]). For 
gas mixtures with L = 1 = a the subcritically bifurcated 
TW solution branch ends on the SOC branch without 
having formed a saddle. Thus, this TW branch remains 
unstable all the way from the bifurcation threshold r osc 
to the end point r* on the SOC solution branch. This 
result is supported by a recent test calculation with a 
267-mode Galerkin expansion. In it the number of 6— 
and (-field modes was large enough to reproduce any 
structural details of the finite-difference solution of the 
field equations. 

The model shows in a small Q-i/'-region — cf the 
shaded areas of Fig. 7 — forwards bifurcating TWs (e.g., 
for Q = 0, L = 1 = a, and k = between the tricrit- 
ical value ip t = —0.3812 and the codimension-two value 
ipCT = —0.2174). However, the r-interval (r osc ,r*) where 
these supercritical TW solutions appear is very small — 
less than 1% of the current heating rate r for Q = — so 
that the initial slope of the bifurcating TW branch is al- 
ways large. On the other hand, in the 267-mode Galerkin 
calculations, done for Q = 0, a = 1, we found only back- 
wards bifurcating TW solutions for L = 1 while runs for 
L = | did show a small ^-range near the CT point with 
supercritical TWs. So the domain boundaries between 
supercritical and subcritical TW states in the L—a—ip— 
Q parameter space are not reproduced quantitatively by 
our 8-mode model. 



2. Stationary convection 

In Fig. 6 we compare the bifurcation properties or our 
model SOC solutions (4.18) with finite-differences sim- 
ulations of the full field equations. We show the Nus- 
selt number N and the mixing parameter M versus r for 
four Dufour numbers (Q = 0, 5, 10, 20) at four different 
separation ratios — note the different abscissa and or- 
dinates scales. Fig. 6 shows that our model reproduces 
the SOC bifurcation diagrams of N and M including the 
Dufour-induced trend towards a supercritical bifurcation 
topology with less convection and less mixing not only 
qualitatively but also semiquantitatively. 

For a more detailed comparison of the SOC field struc- 
ture of the model with that one of the numerically sim- 
ulated solution we show in Fig. 9 vertical profiles of the 
first lateral Fourier modes n = 0, 1, 2. The velocity field 



is well approximated by our first lateral Fourier mode. 
Higher modes not contained in our model carry at most 
5% of the convection, e.g., for the strongly nonlinear state 
at r = 3.5. Of course, their relative contribution in- 
creases with r. At r = 2, slightly above the saddle of 
the SOC-solution, our model reproduces the exact re- 
sults for the temperature field with the same accuracy 
as for the velocity field. Here, increasing the Rayleigh 
number leads to a building up of a plateau in the first 
mode. This structure cannot be reproduced by our model 
since it takes only a cos -kz into account. On the other 
hand, the model approximates the zeroth lateral mode in 
a nearly perfect manner. Higher modes in the tempera- 
ture field contribute with about 10%. The largest errors 
occur in the concentration field. While they are still ac- 
ceptable at r = 2 there are structural differences to be 
seen in the strongly nonlinear state at r = 3.5. There the 
exact profile shows a local maximum structure caused by 
the combination of the terms 1, cos -kz and cos27rz. The 
latter one is absent in our approximation. In contrast 
to this first lateral mode our model reproduces the ze- 
roth lateral mode of the concentration with the accuracy 
obtained for the velocity and temperature field: Apart 
from quantitative errors (of about 15% for r = 3.5) the 
model also reproduces the inversion of the central con- 
centration gradient as the numerics do. For nonlinear 
stationary convection in ip < 0-mixtures we can there- 
fore state that the stable layering of the concentration in 
the purely heat conducting state is inverted in nonlinear 
states. This leads to a strong mixing of the fluid in the 
bulk of the fluid near z = 0. The profiles of our model 
show similar qualities and defects for positive ip. There 
the nonlinear states invert the unstable layering of the 
conductive state into a stable layering. We finally stress 
that, e.g., the third lateral concentration mode that is 
not shown in Fig. 9 for the sake of clarity contributes 
with 20% relative to the first mode. 

Altogehter the model predicts the vertical profiles of 
the first lateral modes reasonably well near the saddle of 
the SOC's. There, even the concentration field is approx- 
imated in an acceptable manner. In strongly nonlinear 
states we can approximate the vertical velocity field and 
that of the temperature with a satisfying quality. 

VI. SUMMARY 

We have investigated the influence of the Dufour cou- 
pling, i. e., the effect that concentration gradients drive 
heat currents or change the temperature field on convec- 
tion in binary fluid layers heated from below. The Du- 
four effect changes the temperature field equation "di- 
agonally" via the term LQip 2 V 2 T that enhances heat 
diffusion and " off-diagnonally" via the term —LQip'V 2 C 
that reflects concentration-induced temperature varia- 
tions. Thus the largest effects are seen for large L and 
large \ ip \. We focussed our main interest on gas mix- 
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tures with Lewis and Prandtl numbers around one. For 
Dufour numbers < Q < 40, i. e., in a range that seems 
to contain physically realistic gas parameters we have 
determined the influence of the Dufour effect, both on 
the linear stability behaviour of the quiescent conductive 
state and on various nonlinear convection properties in a 
range of experimentally relevant Soret couplings ij). To 
that end we have used an eight-mode Galerkin approxi- 
mation that describes convection in the form of straight 
parallel rolls subject to realistic vertical boundary condi- 
tions. Its nonlinear properties were compared with some 
selected numerical simulations of the full hydrodynamic 
field equations and the linear ones were compared with 
numerically or analytically exact stability analyses. 

LINEAR PROPERTIES — The influence of the Du- 
four effect on the stability of the conductive state can be 
summarized — see also [14] — as follows: (1) It destabi- 
lizes (stabilizes) the conductive state against the growth 
of stationary (oscillatory) convection. (2) The critical 
curves r c stat (ij)) and r c osc (ij)) are shifted with increasing 
Q towards more negative ij). (3) Thus the V'-range with 
stationary (oscillatory) instabilities grows (shrinks) with 
Q. (4) The Hopf frequency decreases with growing Q. 
These properties are well reproduced by our few-mode 
model. In agreement with [15] we found that the crit- 
ical stationary wave number k c stat (p) in the exact sta- 
bility analysis is governed by only a single parameter p 
that also determines the shape of the stationary marginal 
stability curve ^R s tab{k; p). Thus, L, Q, and ip enter 
only via the scaling combinations p = l^i + q^2^ 1+ ^ and 
S = (l + Qip 2 )(l + ip) + ip/L. The exact value for which a 
stationary zero wave number instability occurs is found 
to be po = and an expansion around this point is 
presented. 

NONLINEAR PROPERTIES — To determine the in- 
fluence of the Dufour effect on nonlinear convective prop- 
erties we investigated in particular the bifurcation be- 
haviour of the flow intensity, of the convective heat cur- 
rent TV" — 1, and of the convective mixing of the concen- 
tration field M as functions of r for several Q. 

Traveling wave convection — (1) The ^-range in which 
TW solutions — stable or unstable — are present is 
shifted to more negative values when the Dufour cou- 
pling increases. (2) Forwards bifurcating TWs disappear 
already for smaller Q while backwards bifurcating ones 
still exist for larger Q when ij) is sufficiently negative. 
(3) However, for L = a = 1 these subcritically bifurcated 
TW solution branches do not develop a saddle before 
they merge with zero frequency with the SOC solution 
branch. Thus, stable TW states on an upper TW solution 
branch beyond a saddle that can be seen in liquid mix- 
tures for L = O(10~ 2 ) und a = O(10) were not found, 
neither in the model nor in the numerical simulations of 
the full field OBE's. The reason for their absence is the 
large Lewis number L = C(l) of gas mixtures. 

Stationary convection — (1) The Dufour-induced 
destabilization of the conductive state against stationary 



perturbations shifts the onset of SOC to lower values of r. 
(2) The range of subcritical SOC bifurcations in the L- 
ij> parameter plane shrinks upon increasing Q: (3) With 
growing Dufour coupling there is a gradual change from 
a strongly backwards SOC bifurcation — e. g. at large 
negative ij) — via a tricritical one to a forwards one and 
the initial slope of a supercritical bifurcation curve de- 
creases. (4) In addition, the Dufour effect reduces the 
convective intensity and with it the mixing of the con- 
centration field at larger r so that the bifurcation curves 
of N — 1, u)^ al , and M as functions of r become flat- 
ter. (5) Furthermore, structural changes in the convec- 
tive fields resulting from either increasing Q or decreasing 
r are the same. (6) Comparison with numerical simula- 
tions of the full OBE's shows that the Galerkin approx- 
imation provides good results for stationary convection 
in gas mixtures also for SOC states that are well above 
onset not only for the Nusselt number but also for the 
convective mixing M. The reason is that the field ex- 
pansion in trigonometric functions is quite effective since 
with L = 0(1) the problem of concentration boundary 
layers is less severe than in liquids. 
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APPENDIX A: FSI STABILITY ANALYSIS 

Here we consider a FSI-mode truncation to study the 
influence of the Dufour effect on convection in binary gas 
mixtures with an idealized free-slip boundary condition 
for the velocity field but with impermeable plates using 
the ansatz 
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Performing a stability analysis of the conductive state 
in the standard way we reproduce the stationary stabil- 
ity threshold r stat (k) and the critical wave number k c stat 
obtained by Hort et al. [14]. But for the oscillatory insta- 
bility we obtain results that differ for finite Q from eqs. 
(3.17-3.19) of ref [14]. We get 
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with a critical lateral wave number determined by a 
third-order polynomial in k c os 2 
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and a Hopf frequency Wjj(fc) given by 
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There we used the following abbreviations 
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In Fig. 10 we present our FSI results which show in con- 
trast to Fig. 2 of [14] the characteristic Dufour-induced 
shrinking of the ^"region with an oscillatory instability 
of the conductive state. 
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FIG. 1. Dependence of stationary critical properties on the 
scaling variable p (3.6). (a) Reduced critical wave number 
(solid line). The dashed line represents the expansion around 
po = 131/34. (b) Reduced scaled stability threshold. 

FIG. 2. Stability properties of a gas mixture (i = 1, a = 1) 
vs. separation ratio tp for different Dufour numbers Q. The 
reduced stationary (solid line) and oscillatory (dashed line) 
stability thresholds r c stat and r% sr , the corresponding reduced 
wave numbers k c stat and k% sc , and the critical Hopf frequency 
uh result from the NSI model for which R° c = 1728.38 and 
k° c = 3.098. 

FIG. 3. Square of the vertical velocity mode X of station- 
ary convection vs reduced Rayleigh number r for different 
values of the Dufour number Q. Parameters are tp = —0.25, 
L = 1, k = fc", and arbitrary a. The long dashed line 
shows X 2 = r — 1 in a pure fluid (V> = 0). For com- 
parison with experiments one should identify X 2 with the 
flow intensity reduced by the pure fluid value at r = 2, i.e., 
X 2 = w 2 max lw 2 max {ij) = 0,r = 2). 

FIG. 4. Onset behaviour of SOC with k = fc" in the L-ij) 
plane for different Dufour numbers. The bifurcation of flow 
intensity vs Rayleigh number is forwards (backwards) to the 
right (left) of the full thick curve of tricritical bifurcations. 
The threshold r,t a t has moved to r,t a t = oo at the thin line. 
Below this curve the convective solution is disconnected from 
the conduction fixed point. In this parameter regime con- 
vection branches out of the conductive state for heating from 
above, r < 0. 

FIG. 5. Structural properties of SOC states in the x-z 
plane perpendicular to the roll axes. Black implies low field 
values and white high ones. The second row shows the lat- 
eral profiles of w (solid line), 209 (dashed line), and 100c 
(dot-dashed line) at mid height of the fluid layer. Increasing 
Q or decreasing r leads to structurally similar fields — increas- 
ing both appropriately does not cause changes. Parameters 
are tp = —0.25, L = 1, and k = 1. 



FIG. 7. Existence boundaries of TW solutions (stable or 
unstable) of the model for a = 1, k = k° c . TWs exist below 
the thick labelled curves. For k = k c otc the existence range 
is larger as can be inferred from Fig. 2. The shaded (white) 
regions show the range of supercritically (subcritically) bifur- 
cating TWs. 

FIG. 8. Vertical profiles of the linear critical Fourier modes. 
Full lines refer to the exact result obtained analytically (see 
Sec. Ill) or numerically with a Shooting method. Dashed lines 
represent the model. Because of the mirror symmetry of the 
critical modes at the mid plane, z = 0, only one half of the 
layer is shown. Left column contains the moduli at the sta- 
tionary instability for tp = —0.3. Right column shows the 
moduli and phases at the oscillatory instability for tp = —0.4. 
The vertical average of tp w (z) has been assigned to the phase 
angle zero. The normalization is always such that | 9 \= R° c 
at z = 0. Parameters are a = 1 = L and Q = 5. 

FIG. 9. Vertical profiles of the lateral Fourier modes n = 
(circles and full lines ), n = 1 (squares and dashed lines), 
and n = 2 (triangles) for stationary convection. Symbols re- 
fer to finite-difference numerical simulations of the full field 
equations and lines represent our model. In the last two rows 
SC denotes the deviation of the concentration from the global 
mean and the dotted line is the conductive profile of SC. Nu- 
merical profiles of, e.g., the n = 3 concentration mode con- 
tributing with about 20% of the first mode are not presented 
for the sake of clarity. Parameters are tp = —0.25, Q = 0, 
L = a = k = 1. 

FIG. 10. Stability properties of a gas mixture (i = 1, 
<7 = 1) vs separation ratio tp for different Dufour numbers Q. 
The stationary (solid line) and oscillatory (dashed line) stabil- 
ity thresholds r c stat and r% sr , the corresponding reduced wave 
numbers k\ tat and k% sr , and the critical Hopf frequency u> h 
are determined approximately for FSI boundaries for which 
i?" = ^-7r 4 and fc" = 7r/-s/2- This figure is meant s a corrigen- 
dum of Fig. 2 of [14]. 



FIG. 6. Nusselt number N and mixing parameter M of sta- 
tionary convection versus Rayleigh number r resulting from 
our model and from numerical simulations of the full field 
equations with L = a = 1, k = k° c . Within a column the 
Soret coupling tp has the value indicated in the figure. Curves 
for different Dufour numbers Q are identified in the legend. 
Unstable branches (dotted lines) of subcriticallly bifurcating 
SOCs can be seen only in our model, namely for tp = —0.25 
and ij) = —0.5. 
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